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The hematite (a-Fe 2 03) (0001) surface: evidence for domains of distinct chemistry 
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Using spin-density functional theory we investigated various possible structures of the hematite 
(0001) surface. Depending on the ambient oxygen partial pressure, two geometries are found to 
be particularly stable under thermal equilibrium: one being terminated by iron and the other by 
oxygen. Both exhibit huge surface relaxations ( — 57% for the Fe- and — 79% for the O-termination) 
with important consequences for the surface electronic and magnetic properties. With scanning 
tunneling microscopy we observe two different surface terminations coexisting on single crystalline 
a-Fe2C>3 (0001) films, which were prepared in high oxygen pressures. 

PACS numbers:68.35.Bs,61.16.Ch,68.35.Md,71.15.Ap 
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Although metal-oxide surfaces play a crucial role for 
several profitable processes, good quality experimental 
and theoretical studies of their atomic structure and 
electronic properties are scarce. For example, a-Fe2 03 
appears to be the active catalytic material for produc- 
ing styrenejj which was substantiated by recent reac- 
tivity studies performed over single crystalline hematite 
model catalyst Jjlms.tJ Other candidate applications are 
photoelectrodeaJ and non-linear optics materials! Nev- 
ertheless, the surface properties of ev-Fe2 03 are basically 
unknown, and also for other metal oxides an understand- 
ing is developed only badly. The reason is the difficult 
preparation of clean surfaces with defined structures and 
stoichiometries, which, as in the case of hematite, can 
require high oxygen pressures not suitable in standard 
ultrahigh vacuum systems. Furthermore, electron spec- 
troscopy techniques and scanning tunneling microscopy 
(STM) are hampered by the insulating nature of the ma- 
terial. We also note that surface-science techniques of- 
ten do not probe a thermal equilibrium geometry but a 
frozen-in metastable state. Theoretical studies, on the 
other hand, have to deal with 3d electrons, oxygen with 
very localized wave functions, a rather open structure, 
unusual hybridization of wave functions, huge atomic re- 
laxations, big super cells, and magnetism. This renders 
an ab initio study of a-Fe2 03 surfaces a most challenging 
investigation. Some theoretical studies of the geometry 
of a-Fe203 (0001) had been performed using empirical 
(classical) potentialsOfl, and Armelao et aln studied the 
electronic structure employing a cluster approach. In this 
paper we report spin-density functional theory (SDFT) 
calculations for a slab geometry (see Fig. [l])-nWe use 
the generalized gradient approximation (GGA)@ for the 
exchange-correlation functional and the full-potential iiru 
earized augmented plane wave (FP-LAPW) methocoE2l 
to solve the Kohn-Sham equations. The STM study was 
performed on a thin hematite film grown epitaxially onto 
a Pt (111) substrate. 

The identification of thermal equilibrium structures of 
surfaces is a prerequisite for an understanding of the en- 
durance, electronic, magnetic, and chemical properties of 
the material. For a two-component material, as Fe2 03, 
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FIG. 1. Upper half of the slab for the unrelaxed 
03-terminated surface. The cross section of the upper half 
of the supercell is noted by the dashed rectangle. The A indi- 
cate the calculated interlayer relaxations in % relative to an 
unrelaxed surface, i.e., a truncated bulk geometry. 



under realistic conditions the surface will exchange atoms 
with chemical reservoirs. Therefore we will analyze the 
Gibbs free energy with its dependence on the chemical 
potentials of the two components. For metal oxides, it 
is obvious that the O2 partial pressure is the practical 
handle on this dependence. Interestingly, neither exper- 
imental nor theoretical studies have identified the com- 
position and structure of surfaces of Fe203 so far, and 
it is unknown how the surface structure depends on the 
oxygen chemical potential. In this paper we show that 
the surface properties change significantly with O2 pres- 
sure and that under typical oxygen pressure conditions 
two chemically distinct domains are likely to exist. 

The a-Fe2 03 crystallizes in the hexagonal corundum 
structure, where the primitive bulk unit cell contains six 
formula units (30 atoms). Along the [0001] axis this 
structure can be viewed as a stack of Fe and O3 layers: - ■ - 
Fe-Fe-0 3 -Fe-Fe-0 3 - ■ •. The (0001) surfaces stable and 
often occurs on naturally grown crystalstil. It exhibits 
an jtmceconstructed (lxl) surface with a hexagonal unit 
celOEi In the total-energy calculations we considered 
the following (likely) surface terminations with (lxl) pe- 
riodicity: Fe-Fe-0 3 - • •, Fe-0 3 -Fe- • •, 3 -Fe-Fe- • •, 2 - 
Fe-Fe-- • •, Oi-Fe-Fe-- • •, where the latter two surfaces cor- 
respond to an oxygen termination with oxygen vacancies. 
For the Fe-03-Fe-- • ■ terminations we considered all four 
possible geometries of the surface Fe, and find that the 
energetically favorable surface is that with the top Fe site 
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occupied according to the bulk stacking sequence. Fig- 
ure [y shows as an example our slab for the 03-Fe-Fe— • • • 
surface study. We find that two of the five candidates 
considered have a particularly low energy and predict 
that both should be present under typical experimental 
conditions. These are the Fe-03-Fe-- • • and the 03-Fe-Fe- 
• • • structures. The low energy of the Fe-03-Fe-- • • sur- 
face is consistent with pcent findingsJpr other corundum 
type materials (Al^Ei! and ^oJB), but for Fe 2 3 it 
had not been identified so far. The other structure which 
we predict, i.e. the termination by O3, is unexpected in 
general (also for other corundum type crystals) because 
this termination suggests a high surface dipole moment 
and thus an electrostatically unfavorable situation. In 
fact, our calculations reveal that this argument, which 
is based on the understanding of the bulk properties of 
Fe 2 03, is too simple for describing the surface. Indeed, 
both low-energy surface structures are stabilized by a sig- 
nificant surface relaxation which implies that the surface 
should be viewed as a "new material" , not just a trun- 
cated bulk. 

For the FP-LAPW calculations we use a kinetic-energy 
cutoff for the plane-wave basis of E^^. = ^ty- This is 
a rather high value and makes the calculations very in- 
volved. However, because of the huge surface relaxations 
found in the course of this study we had to use rather 
small muffin-tin spheres (Rf^ = 0.95 A , R^ T = 0.74A ) 
and therefore a large value for E^ l ax was mandatory to 
ensure good numerical accuracy. Details of the calcula- 
tions are presented elsewhere. The bulk as well as the sur- 
face calculations were performed with a hexagonal unit 
cell and a uniform k-point mesh with five points in the 
irreducible part of the Brillouin zone. All parameters 
defining the numerical accuracy of the calculations were 
carefully tested. For the tplk we obtain the following 
results (experimental valuesll3 in parentheses): ao=5.025 
(5.035) A, c =13.671 (13.747) A, z(Fe)=0.357 (0.355) A, 
x(O)=0.308 (0.306) A. Our calculations give an antifer- 
romagnetic ordering with a local moment of M = 3.39/iB 
(considering the contribution from the muffin tin only). 
The spins are pointing parallel to [0001], and in the Fe-Fe 
double layers arc aligned parallel. But between neighbor- 
ing double layers (separated by an O3 layer) the spins are 
antiparallel. Our result is similar to recent bulk calcula- 
tions using the augmented spherical wave methodED, but 
the magnetic moment rare-smaller than the experimental 
values M = 4.6-4.9/ib"'"- This differencej-was blamed 
on difficulties of the experimental analysisHEi-, The heat 
of formation per Fe 2 3 unit is 8.024 (8.48)B eV. The 
slabs consist of six O3 layers and ten, twelve, or four- 
teen Fe layers, depending on the surface termination to 
be studied (see Fig. [j] as one example). All atoms of the 
slab are allowed to relax and no symmetry restrictions 
are applied. The antiferromagnetic ordering is found to 
remain up to the surface. 

The Gibbs free energy Q, of the slab at temperature T 
and partial pressure p is given by 
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FIG. 2. Surface energies of different Fe2Os(0001) surface 
terminations, no(gas) is the chemical potential per oxygen 
atom of molecular O2. The allowed range of ^o-fio(gas) is 
indicated by the vertical dotted lines, where the left one cor- 
responds to strongly Fe-rich (i.e., oxygen-poor) conditions, 
and the right one corresponds to strongly oxygen-rich condi- 
tions (i.e., high oxygen gas pressure). Full lines show results 
for relaxed geometries, and dashed lines give for comparison 
results for unrelaxed surfaces. 
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where E total is the total energy of the slab, /z Fe the chem- 
ical potential of iron, and /io the chemical potential of 
oxygen. A Fo and No are the number of iron and oxygen 
atoms of the supercell. For typical pressure and temper- 
ature, the pV and TS terms in Eq. [l] can be neglected. 
The Fe and O chemical potential are not independent, 
because they are related to each other by the existence 
of the Fe 2 03 bulk phase. This gives 



n = # total - ^A FcMFo2 o 3(bulk ) + (^A Fc - A g )mo , (2) 



where /iFe 2 3 (bulk) is the total energy per bulk Fe 2 03 
formula unit (with 2 Fe and 3 O atoms) . 

In Fig. H we show our results of Eq. ||for various (1 x 1) 
geometries, where 7 is the Gibbs free energy per surface 
area. The meaningful range of /lo-IJ-o(gO'S) is limited by 
the conditions that the chemical potential of Fe has to 
be smaller than that of an atom of bulk iron, and the 
chemical potential of oxygen has to be smaller than that 
of an oxygen atom of 2 . Otherwise an iron or oxygen 
condensate will form at the surface. The allowed range 
of the oxygen chemical potential, ^o-fJ-o{g a s), is marked 
in Fig. U by the vertical dotted lines. 

The results show that the Fe and the O3 terminated 
surfaces are particularly stable. The relaxation energy 
is significant, in particular for the O3 termination. In 
fact, from these results we predict that under typical ex- 
perimental conditions the surface should consist of two 
domains, with these two structures. The O3 termina- 
tion has not been discussed before but for the Fe-03- 
Fe-- • • termination semi-empirical studies had been per- 
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TABLE I. Interlayer relaxations at the Fe-03— • • and 
03-Fe- • • surface relative to the corresponding bulk spacings. 
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formed. Our surface energy for the relaxed geometry is 
1.52 J/m 2 = 94.58 meV/A which is in close agreement 
with the semi-empirical studies of of Mackrodt et alu, 
who obtained 1.53 J/m 2 , and Wasserman aln, who ob- 
tained 1.65 J/m 2 . Also other corundum- type materials 
have a comparable surface energy for the (0001) surface: 
For a-A^G^Manassidis et al!E3 obtained 1.76 J/m 2 , and 
Rohr et al&B obtained for Cr2C>3 a value of 1.60 J/m 2 . 

Table | collects the calculated surface relaxations for 
the Fe-terminated surface. We obtain a relaxation pat- 
tern which alternates when going away form the surface. 
For the Fe terminated surface, the comparison with pre- 
vious (serruempirical) studies reveals that a model using 
static ionsu is .Hot appropriate, but when some polariza- 
tion is allowedQ, the system is apparently well described. 
In fact, our calculation shows that at the surface the cova- 
lent contribution to bonding is enhanced, which is mainly 
due to a hybridization of O 2p and Fe 3d orbitals. The 
workfunction of the Fe terminated surface is calculated 
as 4.3 eV (before relaxation it is only 3.1 eY). Recently, 
in a LEED analysis for Cr 2 03 Rohr et al£3 determined 
a 60 % reduction for the first Cr-03 interlayer distance, 
which is in surprisingly close agreement with our result 
for Fe2C>3 (we obtain 57 %). A detailed list with the ge- 
ometry data which we obtained will be given elsewhere. 
Here we only note that the O3 triangles of the first oxy- 
gen layer, i.e. the second layer from the surface, exhibit 
a clock-wise rotation by 2° without breaking the C3 sur- 
face symmetry. No LEED analysis for the Fe-terminated 
(0001) surface of Fe2C>3 exists so far. 

The O3 terminated surface (see Figs. I and |) is in- 
deed a very unexpected system and only stabilized by 
huge and unusual relaxations which are collected in Ta- 
ble [jj The top layer O3 triangles undergo a significant 
rotation (by 10°) without breaking the C3 surface sym- 
metry. The oxygen atoms remain nearly planar. Simi- 
larly to the Fe-terminated surface, a very notable feature 
is the contraction of the first subsurface Fe-Fe interlayer 
spacing, followed by an expansion of the next Fe-03 spac- 
ing. The huge interlayer relaxations and the O3 rota- 
tional reconstruction provide the important contribution 
to the decrease of the surface energy and stabilize the O3- 
terminated surface. As mentioned above, the stability of 



this surface cannot be understood in terms of a simple 
ionic model, in which the surface oxygen atoms where 
negatively charged. Indeed, according to our calculations 
the ionic character of the O atoms is decreased at the O3- 
tcrminated surface, and the covalent interaction between 
O 2p and Fe 3d states is enhanced significantly. This is 
also reflected in a noticeable decrease of the work function 
upon relaxation (from 8.3 eV to 7.6 eV). 

The surface relaxations and change in the nature of 
bonding at the surface goes hand in hand with a dra- 
matic change in the magnetic properties. Here, we dis- 
cuss the integral over the magnetization density consid- 
ering the contribution from inside the muffin-tin spheres. 
Therefore the direct values should not be taken too lit- 
erately. But differences in moments between different 
layers and between the relaxed and unrelaxed situation 
are indeed meaningful. For the Fe-terminated surface we 
find that the magnetic moment of the topmost Fe layer 
is reduced against the bulk value by about 7 % (namely 
AM = -0.24/x B ). 

The electronic and magnetic properties of the O3- 
terminated surface are unprecedented. We find that the 
subsurface iron layers change their character significantly. 
The local magnetic moments of the first and second Fe 
layer are reduced to less than 50 % of the bulk values 
(AM = -1.79 and -1.75 /ib respectively), and again, a 
significant fraction of this reduction only occurs upon 
surface relaxation and reconstruction. Interestingly, also 
the top O3 layer now attains a noticeable magnetic mo- 
ment (M=0.20 hb per atom). In the bulk, the magnetic 
moment of oxygen atoms is zero. Our calculations show 
that the surface states for the 03-terminated surface are 
of Fe 3d character. They are partially occupied by the Fe 
3d electrons formerly with different spin, which results in 
the decrease of the local magnetic moments of the iron 
atoms. Our calculations also show that the surface state 
electrons of the iron atoms reaching through the topmost 
O layer may noticeably affect the surface reactivity. 

The STM experiments were performed in an Jiltra- 
high vacuum system described in detail elsewhereni It 
is equipped with a sample transfer mechanism and a sep- 
arate preparation chamber for performing high pressure 
oxidation treatments. Single crystalline a-Fe2O3(0001) 
films were grown onto a clean Pt(lll) substrate surface. 
In a first step about 10 nm thick Fe304(lll) magnetite 
films were prepared by repeated deposition of iron and 
subsequent oxidation for 2 min at temperatures around 
950 K in 10~ 6 mbar oxygen partial pressures. Then a 
final oxidation at T=1100 K in 10~ 3 mbar oxygen par- 
tial pressure was performed for 10 min. The structure 
formed under these conditions was quenched by cooling 
down the sample to room temperature in the oxygen at- 
mosphere. Then the oxygen was pumped off and the 
sample was transfered back into the analysis chamber 
and immediately studied by STM. These films exhibited 
no Augep-jcontamination signals and sharp (lxl) LEED 
patternsE3. X-ray photoemission spectroscopy revealed 
only Fe 3+ and no Fe 2+ species, indicating the formation 
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FIG. 3. Constant current image of an a-Fe203 (0001) sur- 
face prepared in 10 -3 mbar oxygen partial pressure over an 
area of 100x80 A 2 (U=+1.3 V, J = 1.3 nA). Below a cross 
section along the white line in the image is shown. Two dif- 
ferent surface terminations can be seen. 



calculations predict two kinds of surfaces (thermal equi- 
librium with (lxl) symmetry), which depend on the 
growth conditions. We find a large inwards relaxation of 
the first layer for the Fe-terminated surface and a huge 
contraction of the interlayer spacing between the Fe sub- 
surface layers for the 03-terminated surface. In addition 
to the relaxations along the [0001] direction, the O layers 
of both surfaces have a plane rotational reconstruction. 
To our knowledge, this is the first finding for such kind 
of special relaxations and plane reconstruction of oxygen 
layers at a metal-oxide surface. We also predict an un- 
usual electronic structure of the 03-terminated surface 
with noticeable presence of states from the subsurface 
Fe layer. This also results in a magnetic polarization of 
the oxygen. With STM we confirm the existence of two 
different surface terminations coexisting on single crys- 
talline hematite (0001) films, which have been prepared 
under a high oxygen pressure of 10~ 3 mbar. 



of a single phase hematite filmEl. 

Figure || shows an atomic resolution STM image of 
the a-Fe2C>3(0001) surface. A dark region and a bright 
region can be seen. They both exhibit hexagonal lat- 
tices formed by atomic protrusions with a periodicity 
of 5 A, which corresponds to the interatomic distance 
within the (0001) Fe-planes of the a-Fe2C>3 structure. 
Below a cross section along the white line in the STM 
image is shown, which displays the 5 A periodicity and 
the two regions separated by a step about 1 A high. This 
step height considerably deviates from the distance be- 
tween two equivalent (0001) surface terminations of a- 
Fe2 03, which are separated by a monoatomic step with 
a height of 2.28 A. We always measure monoatomic steps 
or multiples of them between equivalent regions on the 
hematite surface. Therefore the two regions in Fig. || cor- 
respond to two different surface terminations that coexist 
on the a-Fe2O3(0001) surface. This is further substan- 
tiated by their different corrugation amplitudes depicted 
in the cross section plot in Fig. ||, which are 0.1 and 
0.2-0.3 A. Based on the apparent topographic height dif- 
ference measured with the STM, we interpret the bright 
region as the Fe-terminated and the dark region as the 
O-terminated surface. About 20 % of the entire sample 
surface are covered by the bright termination as deter- 
mined by numerous large-area scans. They form islands 
preferentially located near step edges on top of the dark 
termination, with sizes ranging between 20 and 200 A. 
We also confirmed that their formation depends on the 
oxygen partial pressure during the preparation (compose 
the discussion of Fig. || above). Measurements of STM 
spectroscopy, work function etc., and a detailed analysis 
of the dependence on oxygen pressure are in progress. 

In summary, we have presented a detailed ab initio 
FP-LAPW study of the stoichiometry and structural re- 
laxations of the a-Fe2C>3 (hematite) (0001) surface. The 
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